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We demonstrate that the experimental data on the temperature dependence of the resistivity p{T) 
and conductivity (T(r) for high-mobility Si-MOS structures, obtained for a wide range of densities 
(~ 1 : 7) at intermediate temperatures, agree quantitatively with theory of interaction corrections 
in the ballistic regime, Tt > 1. Our comparison does not involve any fitting parameters. 

Recently [1], we compared the resistivity p{T) data 
for high mobility Si MOSFETs with the interaction cor- 
rections calculated to the first order in T [2], and found 
a quantitative agreement over a wide range of densities 
at intermediate temperatures (in the ballistic regime 
Tt > Ti/kB) [1,2]. In the comparison, we used indepen- 
dently determined values of the renormalized g-factor 
and effective mass [3]. 
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In the Comment by Shashkin, Dolgopolov, and 
Kravchenko (hereafter SDK) [4], several claims have been 
made; we respond to the claims below: 

(1) SDK claim that our paper [1] contradicts the ear- 
lier publication [5]. In fact, Ref. [1] focuses on electron 
transport in the ballistic regime {Tt > h/ks), whereas 
Rcf. [5] considers the single-particle interference in the 
diffusive regime {t^p ^ r); it has been shown that the 
weak localization and interaction corrections in the dif- 
fusive regime are irrelevant to the "metallic" p{T). This 
conclusion does not contradict our analysis in Ref. [1] . 

(2) SDK are concerned about our comparison with the 
theory of interaction corrections by Zala, Narozhny, and 
Aleiner (hereafter ZNA) [2] in terms of Ap{T), and claim 
that "there is no linear- in T interval" in our cr(T) data. 

Wherever the T-dependent corrections to the conduc- 
tivity are small (Acr(r)/cr <c 1), one can treat Ap(T) 
and Aa{T) on equal footing. At sufficiently large densi- 
ties {n > 3.5 X 10^^ cm^^), this condition is fulfilled over a 
wide interval of temperatures {h/2TrT < Tiow, Thigh ^ 
Tp, Thigh — Tiow ~ a few K); this enables to detect un- 
ambiguously the linear- in-T regime of a-{T). In Ref. [1] 
we analyzed resistivities, whereas here we compare the 
data for both samples from Ref. [1] with the theory in 
terms of cr(T). Figures 1, 2 show that the cr{T) data are 
in a quantitative agreement with the interaction correc- 
tions to the conductivity a = aD + A.a{T) (solid lines). In 
the calculations, we have used the values for the density- 
dependent renormalized effective mass m*/mb and g- 
factor g*/gb = (x*/x^)("*b/™*)7 determined in indepen- 
dent experiments [3] ; the value of r for each curve is ob- 
tained from the experimental ao = {,p{T — > 0))^^ data 
via the extrapolation procedure described in Refs. [1,2]. 
Once the r value (and, hence, the Drude conductivity 
ctd) is fixed, the slope da/dT is obtained with no free 
parameters. 
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FIG. 1. Comparison of the (t(T) data (symbols) for 
sample Si22 reproduced from Ref. [1] with calculated 
quantum corrections [2] (solid lines). Densities, from top 
to bottom, are (in units of lO^^cm"^): 21.3, 18.9, 16.5, 
14.1, 11.7, 10.5, 8.1, 7.5, 6.9, 6.3, 5.7. 



In the comparison, in accord with Ref. [6] , we consid- 
ered only the temperature range h/2TTT < T ^ Tp, and 
ignored the lower-T and higher-T data. As was discussed 
in Rcf. [1], the temperature range where the agreement 
holds, shrinks towards low densities. This might be ex- 
pected because, as the density is decreased, the p{T) de- 
pendence becomes so steep that the temperature changes 
of Ap/ po remain small only within a narrow range of T. 
The agreement of the data with the interaction correc- 
tions in Figs. 1, 2 holds over a reasonable range of tem- 
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peratures and in a wide range of densities (1:7), but not 



too close to the critical density value ric 
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G{2kF)] + Ikp/qs), and 0(2^^) = 0.2236. By re-writing 
T in units x = Tt/Tl [2], and a in units of /{-kTl)^ we 
obtain 
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FIG. 2. Comparison of the <y{T) data (symbols) for 
sample Sil5 [1] with the interaction corrections [2] (solid 
lines) and the "screening" theory [7] (dashed lines) . Den- 
sities are indicated in units of lO^^cm"^. 

Once again, no matter how good (or bad) the agree- 
ment at low densities is (this regime is beyond the appli- 
cability of the theory [2] anyway), the main conclusion 
of our paper [1] remains the same: the "metallic" behav- 
ior of p(T) in high-mobility Si MOSFETs in the regime 
p <C /i/e^ is in quantitative agreement with the predic- 
tions [2] based on considering the interaction effects in 
the ballistic regime. 

(3) SDK claim that in the range of densities n > 
2 X lO^^cm"^, the p{T) changes can be described by the 
"traditional screening" theory by one of the authors [7] . 

The "screening" theory [7] predicts the dependence 



<5ct(T)/ct(0) = 1 - 2C{a)Cin){T/Ef 



(1) 



with C{a) = 21n2 (or (8/3)ln2) for charged impurities, 
a = —1 (or for surface roughness, a = 0) scattering, 
respectively, C{n) = F{2kF){l - G{2Kf) /{F(2kF)[l - 



da/dx « -2C(a,n) = -4C(a)C(n). 



This expression differs from that derived in Ref. [2]: 



(2) 



da/dx « 15 



1 - -t(x) 



1 - 0/(2;) 



(3) 



where the functions f{x) and t{x) are defined in Ref. [2]. 

The difference between the predictions of the "screen- 
ing" theory [7] and the theory of interaction corrections 
[2] is illustrated in Figs. 2,3. Note that the difference in 
the slopes da/dT in Fig. 2 is a factor of 2 at the lowest 
density n — 2.2 x lO^-'^cm"^ (to calculate the interaction 
corrections to a, we have used the experimental values 
of Fq^u) and m*(n) measured by us in independent ex- 
periments [3]). The theoretical predictions for the slopes 
da/dx as functions of the density are shown in Fig. 3. 
The two density dependences cross each other at differ- 
ent n depending on the specifics of impurity or interface 
scattering in the screening theory. Adjusting the latter 
parameters does not improve agreement between the the- 
ory [7] and our experimental data, if a sufhciently wide 
range of densities is considered. At the same time, the 
ZNA theory [2] is consistent with the data in the wide 
density range without any adjustable parameters. The- 
oretical aspects of the difference between approaches [7] 
and [2] have been clarified in Ref. [2]. 
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FIG. 3. Comparison of the slopes da/dx, calculated 
from Ref. [7] for two models of scatterers (a = — 1 and 
Q = 0), and from Ref. [2] with Fq and m* values deter- 
mined in [3]. In the latter case, the slope slightly varies 
with the temperature (compare x = 1 and x = A), due 
to the functions fix) and t{x) in Eq. (2). Symbols depict 
the slope da/dx of the data shown in Figs. 1, 2. 
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(4) SDK suggest that the low-temperature deviation of 
our p{T) data from the theory, observed for several sam- 
ples, was caused by electron overheating. Apparently, 
they overlooked the corresponding discussion in our pa- 
per (p. 3, left column [1]), which analyzed and rejected 
this possibility. For the range of relatively high densities 
n ^ TT-c wc explored in Ref. [1], the deviations from the 
theory occur at high temperatures, where overheating is 
highly unlikely. Moreover, the low temperature part of 
the data at T < %/2'kt is irrelevant to the linear-in-T 
dependence [6]. 

(5) SDK claim that the range of densities explored in 
[1] is "not related to the anomalous increase of p{T) with 
temperature". The term "metallic" refers to the sign 
of dp/dT, not to its magnitude; however, the "anoma- 
lous metallic" behavior is more commonly referred to the 
strong changes of p{T) [6]. We emphasize that, though 
we have limited our quantitative analysis in [1] to the 
small changes Ap{T)/p <g; 1, where the theory [2] is ap- 
plicable, the overall increase of the resistance with T is 
large in the whole explored range of densities: p increases 
with temperature by a factor of 2 even at n as high as 
15 X 10"cm-2 [5]. 

Note that this claim by SDK is inconsistent with their 
own publications: indeed, in Ref. [8], they fitted the lin- 
ear (t(T) dependences at low densities and low temper- 
atures (T < IK) with the theory [2], whereas on pages 
2,3 of Ref. [9], one of the authors concluded on the basis 
of analysis of cr(T) over the same ranges of densities and 
temperatures that "the linear p{T) dependence is differ- 
ent from the dependence seen at higher temperatures 
and higher densities" . 

(6) SDK claim that the theory [2] "is inapplicable at 
high density n > 5 x lO^^cm^^, where the surface rough- 
ness scattering becomes dominant". In fact, the theory 
[2] assumes the short-range scattering (for more details, 
see Ref. [10]). This requirement is well justified for our 
samples over the entire range of n = (2.2 — 30) x 10^^ 
cm~^ as it has been shown in Ref. [1]: the ratio of the 
"large-angle" scattering time t to the quantum life time 
(or "all-angle" scattering time) Tq = h/{2nTj:)) is close 
to unity (here. To is the Dingle temperature). Besides, 
it is expected that the difference between the effects of 
short-range and long-range disorder on p(T) is negligible 
atT<0.1TF [6]. 

(7) Criticizing our analysis of p{T) at "high" densities 
n ric, SDK give preference to comparison with the 
theory at low densities n ric and low temperatures. 
An attempt to analyze p{T) in this regime on the basis 
of the theory [2] , undertaken by SDK in Ref. [8] , suffers 
from several drawbacks: 

(i)At low densities n w ric and low temperatures, such 
as explored in [8], T is of the order of h/2'KT, and the 
requirement for the ballistic regime is not satisfied (see, 



e.g. [6]). In this case, one cannot neglect the diffusive 
terms and crossover functions, as it was done in Ref. [8]. 

(ii) In the vicinity of ric, the temperature dependence of 
<j{T) is strongly affected by the procedure of sample cool- 
ing [11], this regime is beyond the scope of theory [2]. 

(iii) In the vicinity of n = ric, the ratio Tg/r = 3, as ob- 
served by SDK [8,12] is clearly anomalous; this indicates 
that the standard models of disorder fail. 

(iv) For low densities n lO^^cm"^ (and, correspond- 
ingly, low values of Tp ^ 2.5K), the "metallic" behavior 
corresponds to the temperature range where the valley 
splitting cannot be ignored. (Note that in our anal- 
ysis [1], the valley splitting A^, -C T and, therefore, is 
insignificant). The sample-dependent parameter A„ has 
not been measured for the SDK's sample, and SDK made 
a doubtful assumption that the valley splitting ^ liiT 
[8] , which contradicts their own Shubnikov-de Haas data. 
Correspondingly, in fitting the above Eq. (3), the number 
of triplet components was incorrectly reduced from 15 to 
7. 

For the aforementioned reasons, the interaction correc- 
tions theory [2] can be hardly applied to the (j(T)-data 
analysis at the densities n « Uc. The "accurate mea- 
surements in the best samples" in Ref. [8] (though with 
mobility 2.9 x lO'^cm^/Vs [8] lower than 4 x lO'^cm^/Vs 
as in Fig. 2) have not been analyzed accurately enough to 
provide a reliable quantitative information on the renor- 
malized effective mass. As the result of several unjustified 
simplifications, the authors of Ref. [8] overestimated the 
increase of the effective mass with decreasing density. 

To summarize, 

• Our p{T) measurements performed over wide 
ranges of densities and temperatures, being com- 
bined with the independent measurements of renor- 
malized Fg{n) and m*{n) [3], enable a rigorous 
test of the theory [2]. The theoretical predictions 
are in a quantitative agreement with the experi- 
mental data, in terms of either p{T) or a{T), over 
the density range where the theory can be applied 
(p < /i/e^, Ap{T)/p < 1), and in the ballistic T- 
range. 

• The alternative attempt of such comparison [8] in 
the regime of low densities n ^ Uc and low tem- 
peratures suffers from several drawbacks, and does 
not provide a reliable quantitative information on 
renormalization of the effective mass. 

• We point out that the predictions of the "screen- 
ing" theory of Ref. [7] differ significantly from those 
of the theory of interaction corrections [2] over a 
wide range of densities, which is of primary impor- 
tance for the problem of "metallic" resistivity of 
high-mobility Si MOSFETs. 
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